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We investigate the properties of satellite galaxies formed in N-body/SPH simula- 
tions of galaxy formation in the ACDM cosmology. The simulations include the main 
physical effects thought to be important in galaxy formation and, in several cases, 
produce realistic spiral discs. In total, a sample of 9 galaxies of luminosity comparable 
to the Milky Way was obtained. At magnitudes brighter than the resolution limit, 
My = —12, the luminosity function of the satellite galaxies in the simulations is in 
excellent agreement with data for the Local Group. The radial number density pro- 
file of the model satellites, as well as their gas fractions also match observations very 
well. In agreement with previous N-body studies, we find that the satellites tend to 
be distributed in highly flattened configurations whose major axis is aligned with the 
major axis of the (generally triaxial) dark halo. In 2 out of 3 systems with sufficiently 
large satellite populations, the satellite system is nearly perpendicular to the plane 
of the galactic disc, a configuration analogous to that observed in the Milk Way. The 
discs themselves are perpendicular to the minor axis of their host halos in the inner 
parts, and the correlation between the orientation of the galaxy and the shape of the 
halo persists even out to the virial radius. However, in one case the disc's minor axis 
ends up, at the virial radius, perpendicular to the minor axis of the halo. The angular 
momenta of the galaxies and their host halo tend to be well aligned. 



1 INTRODUCTION 



A successful theory of galaxy formation needs to explain a 
wide variety of physical phenomena across many scales. On 
the Galactic mass scale, the problem of how small satellites 
form around galaxies such as the Milky Way poses two inter- 
esting, and as yet not fully answered questions. These are: 

1) What determines the number of Milky Way satellites 
and the shape of their luminosity function? 

2) Why are the Milky Way's satellites aligned on a great 
circle in the sky and why is this great circle nearly perfectly 
perpendicular to the disc of the Milky Way? 

At present, it is unclear whether the luminosity function 
of the satellites of the Milky Way and their peculiar spatial 
distribution are unique to the Galaxy, or whether these are 
generic features that arise during the formation of a typical 
galaxy. By virtue of their small size, extragalactic satellites 
are difficult to observe and most studies are hampered by 
either flux limits, completeness, or sample size. However, re- 
cent advances using large sky surveys have resulted large 
samples - albeit with only very few satellites per primary - 
and have shed some light on these questions. 

In a universe in which the gravitationally dominant 
component is cold dark matter (CDM), the existence of 



small-scale power in the initial conditions causes the early 
collapse of matter; cosmological simulations show that small, 
dense, CDM haloes are able to form at early times. These 
low mass haloes grow by successive mergers and smooth ac- 
cretion g iving rise to the large-scale structures we see to- 
day (e.e iFrenket alJ Il985t iNavarro. Frenk fc White! 1199^ 
IWechsler et alJl2004t Izhao et al]l2003l) . High resolution N- 
body simulations of the growth of CDM halos have shown 
that the dense cores of merging clumps often survive the 
disruptive effects associated with mergers and remain as dis- 
tinct substructures (or subhalos) embedded within a larger , 
smooth main halo llKlvpin et alJll999t iMoore et al.lll999T) . 
Although the total mass attac hed to subhalos is only of order 
~ 10% of the total halo mass jGhigna et al.ll99StlGao et alJ 
12004 ISprineel et alJl20oitlStoehr et alJl2002l) . both N-body 
simulations and semi-analytical calculations of the assembly 
of dark matter haloes based on the extended Press-Schechte r 
theory llBond et alJll99ltlBowe3ll99lllLacev fc Colelll993f) . 
show that many more embedded substructures survive than 
there are visible satellites in galactic halos (Kauffmann, 
Whit e & Guiderdoni 1993; IMoore et al1ll999l : iKlypin et ail 
1999). 

Many authors have argued that this so-called "miss- 
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ing satellite problem" poses a severe challenge to the CDM 
cosmology and e ncourages the study of alternative forms 
of dark matter dCraie fc Davis! hoOll: iMoore et alJ 120001: 
ISpereel fc SteinhardtfeOOd lYoshida, et al.ll2000f) or of differ- 
ent c osmological initial conditions (iKam ionkowski & Liddlc 
2000). Other authors have argued that the paucity of satel- 
lites in the Local Group is a natural byproduct of the 
physics of galax y formation that r egulate the cooling of gas 
in small halos ( Kauffmaim^^lJ^^dj Bullock, Kravstov 
& Weinberg 2000: iBenson et al.ll2002al) or merely the re- 
sult of a misidentification of substr uctures in the simul ations 
with satellites in the Milky Way ijStoehr et aljl2002>). The 
semi-an alytical model of galaxy formation of IBenson et al] 
(2002b) included a detailed treatment of the reionization of 
hydrogen in the early universe which, by altering the ther- 
modynamic state of primordial gas, inhibits the formation 
of small satellite galaxies. In this model, the increase in the 
entropy of the gas has two effects: it inhibits the cooling of 
new gas into small halos and it delays star formation in gas 
that has already cooled. The end result is a satellite popu- 
lation with about the observed numbers seen in the Local 
Group and with a luminosity function that matches the faint 
end of the observed function but not its bright end where 
the model fails to produce enough large, LMC-like satellites. 

The spatial distribution of satellite galaxies poses 
another interesting problem within the CDM paradigm. 
Whereas N-body simulations show that the substructures 
that survive within halos tend to be nearly spherically dis- 
tributed, the galactic satellites of the Milky Way are con- 
fined to a highly flattened structure , a puzzling fact first rec- 



ogniz ed 30 yea rs ago jLvnden-Bel]|il976l. iKunkel fc Demersl 
1979: see also lLvnden-BerJll982h . Kroupa, Thies & Boily 
(2005) drew attention to this discrepancy and concluded 
that the anisotropic alignment of the Galaxy's satellites con- 
tradicts the CDM model. Recent work, however, has shown 
that the satellite galaxies do not populate a random selec- 
tion of subhalos but are preferentially found in a biased 
subset which is arranged in a flattened configuration. This 
bias partly reflects the preferential infall of the most mas- 
sive dark matter clumps along the filaments of the cosmic 
web. This phenomenon is clearly seen ftp various degrees ) 
in the N-body simulat ions analyzed bv |Kanget alJ (f-2005) , 
iLibeskind et alJ ll2005l) and lZentner et all ll2005l) (and in the 
cluster mass simulations of lKnebe et aljl2004l) . These stud- 
ies differ in the precise way in which subhalos are identi- 
fied with satellites but they all agree that flattened satellite 
configurations such as that seen in the Milky Way are not 
uncommon. In particular, Libeskind et al. and Zentner et 
al. followed the formation of satellites by applying a semi- 
analytic galaxy formation model to high-resolution N-body 
simulations. Both studies found not only fattenings consis- 
tent with that seen in the Milky Way, but also that the long 
axis of the flattened satellite distribution tends to be aligned 
with the long axis of the parent dark matter halo. 

Beyond the Local Group, a number of studies have 
claimed correlations between the orientation of central 
galaxi es and the distribution of their satellites. iHolmberal 
( 1969) first identified a lack of satellites in the plane of a 
coadded sample of c entral galaxies out to a projected ra- 
dius of r p < 50 kpc. IZaritskv et al.l (Il997t) found evidence 
for this "Holmberg effect" but only on much larger scales 
(300 < r p < 500 kpc). The reality of the Holmberg ef- 



fect remains controversial. Early authors claimed that if an 
anisotr opy exists at all, it is, even if significant, at best 
small (lHawlev fc Peebles! Il975l ISharp. Lin fc White! Il979l 
iMacGillivrav et alll982tL An enhancement perpendicular to 
the disc was recently inferred in the distribution of sa tellites 
identified in the 2dFGRS bv lSales fc Lambasl ll2004l) but a 
private communication quoted in lYang et al] J2005I) indi- 
cates that the original analysis was incorrect and that the 
enha ncement is , in fa ct, a long the disc, not perpendicular 
to it. iBrainerdl <2005l) and lYang et al] ]2005h have claimed 
to see an alignment of satellites in the SDSS 1 in the op- 
posite direction to Holmberg's, that is along t he plane of 
the ga laxy disc rather than orthogonal to it. lYang et al] 
(2005) found that the planar distribution is detectable only 
in red satellites but the blue population is consistent with an 
isotropic distribution. These apparently conflicting observa- 
tional studies all seem to suggest anisotropic distributions of 
satellites but they disagree on how the overall distribution 
of satellites is oriented relative to the central galaxy. 

In principle, simulations are an ideal way to investigate 
this sort of issues. Since th eir introductio n to cosmology in 
the 1970s and early 1980s <lPeeblesHl97ll: Aarseth, Turner, 
& Gott 1979, Frenk, White & Davis 1983), N-body simula- 
tions have been extremely useful in revealing how cosmic 
structures emerge out of small primordial perturbations (see 
Springel, Frenk & White 2006 for a review). To investigate 
questions such as the alignment between satellites and 
central galaxies, however, it is necessary to follow not only 
the evolution of dark matter, but the coupled evolution of 
the baryonic component as well. Until recently, progress in 
this area was hampered by the inability of hydrodynamic 
simulations to produce realistic discs from CDM initial 
conditions. Without some form of feedback to prevent most 
of the gas from cooling into subgalactic fragments, the out- 
ward transfer of orbital angular momentum to the dark halo 
as these fragments merge r esults in discs that are much too 
small (iNavarro fc Benzl Il99lt IWeil. Eke fc Efstathioul 
| l99St ISommer-Larsen. Gelato fc Vedell ~ Il999l : 

lEke. Efstathiou fc Wrigtrdl2000D . 

The "disc angular momentum problem" has recently 
been overcome, at least partially, in a number of sim- 
ulations which include plausible forms of feedback and 
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simulations carried out by 'Oka moto et al] (120051) . Specifi- 
cally, we investigate the properties of satellite galaxies or- 
biting central galaxies of mass similar to that of the Milky 
Way. We derive the satellite luminosity function over a lim- 
ited, but still interesting, range of luminosity. We search for 
anisotropy in the satellite galaxy distribution, and study the 
alignment of satellite systems with their central disc, as well 
as the alignment of the disc with the its host dark matter 
halo. 

This paper is organized as follows. In Section 2 we de- 
scribe the simulations we have used, as well as our method 
for selecting complete satellite samples. In that section we 
also derive the satellite galaxy luminosity function and in- 
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vestigate the gas fraction of the largest satellites. We present 
our analysis of relative shape alignments in Section 3 and of 
angular momentum alignments in Section 4, and conclude 
in Section 5. 



2 IDENTIFYING GALAXIES AND 
SATELLITES 

In this section, we briefly describe the simulations that we 
have analyzed and the methods that we have developed in 
order to identify central and satellite galaxies. 

2.1 The simulations 

We have analyzed two simulations of galaxy formation, both 
carried out usin g the parallel PM- TreeSPH code G ADGET 2 
jSnringel l2005h . as modified by lOkamoto et"aH l)200it . 
GADGET2 calculates the evolution of dark matter using N- 
body techniques and the evolution of gas using smooth par- 
ticle hydrodynamics (SPH). 

Initially, the two simulations followed the evolution of 
dark matter in a cosmological cubical volume of length 
= 35.325/i _1 Mpc in a ACDM model with cosmological 
parameters A=0.7, f2 m =0.3, Ho = 70 km s _1 and <Tg=0.9. 
For the first simulation (hereafter SD), a region around a 
suitably chosen dark matter halo was identified at the fi- 
nal time and the simulation was run again, this time adding 
many more dark matter particles, as well as SPH particles, in 
the region of the halo, and perturbing these with additional 
high frequency power drawn from the same power spectrum 
of the ori ginal simulation following the general method out- 
lined bv lFrenk et, al.l l|l996f) . Since the goal of this simulation 
was to investigate the formation of a galactic disc, the halo 
chosen for resimulation was selected to have a quiet recent 
merger history, with no major mergers since z ~ 1. A pre- 
liminary semi-analytical calculation ap plied to the merge r 
tree of this halo, using the methods of iHellv et alJ (2003), 
indicated that a disc galaxy was likely to form in this halo. 
The high resolution region enclosed a spherical volume of 
radius 0.9 Mpc around the halo at z = 0. 

Simulation SD did indeed form a re asonably realistic 
disc, as discussed by Ok amoto et al (in their "shock- 

burst" model). Encouraged by this success, we ran a second 
simulation (hereafter SR) with the same code, this time pop- 
ulating several regions of the same volume with high reso- 
lution dark matter and gas particles. The high resolution 
regions consisted of a large sphere of radius 5/i~ 1 Mpc and 
four smaller overlapping spheres each of radius r=l/i _ Mpc. 
This arrangement ensured coverage of all the large galaxies 
that formed out to the virial radii of their halos. The same 
cosmological parameters were used in both simulations, ex- 
cept that the baryon density was taken to be J7b=0.040 in 
SD, and slightly larger, Ob=0.044, in SR. In both simula- 
tions, the mass per particle was ~ 2.6 x 10 6 for gas and 
~ 1.7 x 10 7 /i _1 Mq for dark matter. The gravitational soft- 
ening was 0.5 and lh" 1 kpc for SPH and high resolution 
dark matter particles respectively, while the minimum SPH 
smoothing length was 0.5 h^ 1 kpc for simulation SD and 
0.39 h~ 1 kpc for simulation SR. 

The various physical proc esses included in our simula- 
tions are described in detail in lOkamoto et all i2005T) . Here, 



we summarize the salient points. The interstella r med ium 
(ISM) is modelled, following lSpringel fc Hernciuisd J2003h . as 
a two phase medium composed of hot ambient gas and cold 
gas clouds in pressure equilibrium. Gas heating and cooling 
rates are computed assuming collisional ionization equilib- 
rium in the presence of a uniform and evolving ultra-violet 
(UV) background, which is assumed to b e generated by hot 
OB s tars and is switched on at z — 6 (Haar dt fc MadavJ 
1996). The cooling rates, which depend on the metallic- 
ity of the gas, are c ompu ted from the tables given by 
ISutherland fc Dopital dl993F) : molecular cooling and other 
forms of cooling below T » 10 4 K are ignored. 

Star formation can occur in a "quiescent" and a "burst" 
mode. In the quiescent mode, gas particles that meet a spec- 
ified density criteria, are turned into stars according to a 
pre-determined probability. Thes e stars form w ith a stan- 
dard initial mass function (IMF; Salpeter 1955). Bursts of 
star formation are triggered by major mergers which are 
identified by tracking large changes in the entropy of the 
gas. In a burst, stars form on a shorter timescale than in the 
quiescent mode and with a top-heavy IMF. This model of 
star formation is m otivated by the semi-analytical work of 
iBaugh et alJ |2005) who argue that only a top-heavy IMF 
in bursts can explain the number density of sub-millimetre 
and Ly man-break galaxies at high redshift. lNagashima et all 
(200 5allbT) argue, similarly, that this model is also required 
to explain the metallicity of the intracluster medium and of 
elliptical galaxies. 

In the Oka moto et alJ i2005h model, the evolution of the 
stellar populations that form is followed in detail, tracking 
both type-II and type-la supernovae. This requires aban- 
doning th e instantaneous recycling app roximation (IRA) as- 
sumed bv lSpringel fc HernauistH2003l) . whereby star forma- 
tion, cold gas cloud formation by thermal instability, the 
evaporation of gas clouds and the heating of ambient gas by 
supernovae explosions all occur simultaneously. Instead, fol- 
lowing each star formation event, supernovae energy and 
metals are injected back into the ISM on a timescale 
which is computed from the mass-dependent stella r lifetime 
iPortinari. Chiosi. fc Bre ssan 1998; Marigo 200l{) and the 
assumed IMF. 

A top heavy IMF in bursts increases the number of su- 
pernovae that explode per unit of mass turned into stars, 
thereby generating stronger feedback. This, in turn, inhibits 
the early collapse of cold gas clouds in small subgalactic ha- 
los, helping to maintain an abundant reservoir of hot halo 
gas. Following the last major merger and accompanying star- 
burst, radiative cooling of hot gas from the reservoir flows 
inwards and settles into a centrifugally supported disc which 
becomes unstable to star formation. This is the key to the 
formation of a realis tic disc galaxy in the simulations of 
lOkamoto etUI (I20051) . 

2.2 Identifying central and satellite galaxies 

We begin the process of finding central and satellite galax- 
ies by identifying "friends-of-friends" (FOF) groups in the 
dark matter, linking together particles whose separation is 
less than 0.2 times the mean interparticle separation, cor- 
responding roughl y to particles within the virialised halo 
iDavis et al.lll985T> . We then identify bound substructures 
("subhalos") in the simulation with the algorithm SUBFIND 
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JSprineel et al.ll200ll) . Using particle positions and veloci- 
ties, SUBFIND calculates the binding energy of each FOF 
group (stripping unbound particles of group membership) 
and then identifies self-bound substructures inside the par- 
ent halo. We consider only halos with ten or more dark mat- 
ter particles, corresponding to a subhalo mass resolution of 
M res ,sub = 1.7 x 10 s ft -1 Mq. For each subhalo, we calculate 
the centre of mass, as well the size, r au b, defined as the rms 
distance of its particles from the centre. 

Substructures within the parent halo are thus identified 
as locally overdense, self-bound, regions in the dark matter 
density field that fall within the high resolution region of 
the simulation. In order to identify clumps of star particles 
as individual galaxies, we associate each star particle with a 
unique dark matter substructure. For each gas and star par- 
ticle, we find the substructure whose centre of mass is clos- 
est and assign the SPH particle to that substructure if it is 
within r su b. Since star clumps tend to be dense and centrally 
concentrated, our results are robust to the choice of subhalo 
radius (for example, increasing r su b by 50% has a negligi- 
ble effect on our results). We also impose a ten star particle 
lower limit, corresponding to M res ,gai ~ 1 X 1O 7 /i -1 M0, on 
the stellar mass of any galaxy. With this resolution limit, 
our galaxy samples become incomplete at magnitudes fainter 
than M v ~ -12. 

We identify "central" galaxies of luminosity similar to 
the Milky Way with galaxies brighter then -20.4 in the V 
band. We calculate their virial radius by growing concentric 
spheres around the galaxy and noting where the mean inter- 
nal density first falls below 200p C rit, where p cr it is the critical 
density. Galaxies that fall within this radius are considered 
satellite galaxies. In order to obtain isolated systems simi- 
lar to the Milky Way, we ensure that no two central galaxy 
candidates are within each other's virial radius. 

Fig. shows an illustrative example of how our algo- 
rithm selects central and satellite galaxies. We plot all the 
star particles within the virial radius projected onto the 
principal planes of the central galaxy, i.e. the planes on 
which the galaxy is face-on and edge-on. In red, we plot 
the star particles that are associated with satellite galaxies. 
There are different numbers of satellites in each projection 
because we are plotting all the star particles in projection 
through the virial radius, which is larger than the plotted 
box. In the bottom right panel, we show the face-on pro- 
jection of the system out to the virial radius. The central 
galaxy has a very extended stellar halo (black points in the 
figure) made up of stars that lie within the r su b of the cen- 
tral galaxy' subhalo and not within the r su b of any satellite 
subhalo. This extended halo makes up only a small fraction 
of the stellar mass of the central galaxy. Below, when we 
calculate the shape of the central galaxy we remove the ex- 
tended halo by considering only the innermost 98% of the 
stars. 



3 THE LUMINOSITY FUNCTION OF 
SATELLITE GALAXIES 

Our simulations produced a sample of 9 central galaxies with 
Mv < —20.4, 8 from simulation SR and one from simula- 
tion SD. The properties of these 9 systems are displayed in 
Table 1. 
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Figure 1. A central galaxy and its satellite system. The black 
dots denote all the star particles associated with the galaxy, while 
the red dots denote the star particles associated with satellites 
located within the virial radius of the central galaxy. The top two 
and the bottom left panels show a zoomed-in projection of the 
distribution of stars near the central galaxy, in the three principal 
planes, i. e. in the planes in which the galaxy is face-on and edge- 
on. The bottom right panel shows a face-on projection out to 
250 kpc. The virial radius, denned by p(< r v i r ) = 200p cr it, is 
marked by the green circle. 



The V-band luminosity function of the satellites in this 
sample is shown in Fig. The luminosity function is nearly 
flat in the range —13.5 > My > —16.5 and drops off sharply 
at fainter and brighter magnitudes. The nominal resolu- 
tion limit of our simulations corresponds to a satellite ab- 
solute magnitude My ~ — 12 and it is important to check 
whether the decline in the luminosity function in the range 
— 12 > Mv > —13.5 is the result of feedback processes, such 
as photoionization or supernova feedback, which affect the 
number of faint galaxies in the simulations, or whether it is 
due to limited resolution. 

To investigate the likely effects of resolution in our es- 
timate of the satellite luminosity function we have used 
the se mi-a nalytic model GALFOR M described by ICole et all 
fcOOdi and lBenson et all i2002aT> . The semi-analytic model 
includes all the standard physical effects present in the simu- 
lation: gas cooling according to a metallicity-dependent cool- 
ing function, star formation, feedback due to supernovae ex- 
plosions, etc. Photoionization is included in an approximate 
way by assuming that gas cannot cool in halos with circular 
velocity v c - 1IC < 6 km s -1 after the as sumed epoch of recom- 
bination, z < 6. iBaueh et alJ i2005t) have shown that this 
simple approximation gives an excellent match to a detailed 
calculation based on the co ncept of a filtering mass llGnedin| 
12001 iBenson et alJl2002bl) . 

We proceed as follows. In the semi-analytic model, it is 
possible to specify the mass of the smallest halo to be con- 
sidered and this allows us to model the effects of resolution 
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Table 1. Properties of the 9 galaxy haloes (gh). We show the virial radius (r v i r ), the dark halo mass internal to this radius (Md m ), 
the stellar mass of the galaxy (iW*), the gas mass of the galaxy (M gas ), the V band magnitude (My), the number of satellites orbiting 
within the virial radius (AT sa t), and the bulge (Mtuige) and disc (AT<i} sc ) mass of the galaxy. 
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in the SPH simulation (iHellv et alJl200St . We recall that in 
the simulation itself we only considered halos with at least 
ten particles, corresponding to M rcSjSU b = 1.7 x 10 s /i -1 Mq. 
We ran a series of semi-analytic models with resolution 
varying by 4 orders of magnitude, from ~ 10 _2 M res>su b to 
10 2 M rcs , su b. In all cases, we found that the shape of the lu- 
minosity function brightwards of Mv = — 12 was essentially 
unaffected. If the resolution is degraded further, then the 
faint end of the luminosity function becomes truncated at 
increasingly bright magnitudes. We conclude from this test 
that our satellite sample is likely to be complete for mag- 
nitudes brighter Mv — — 12 and that our estimate of the 
satellite luminosity function in this regime is unlikely to be 
affected by resolution. 

In Fig. |21 we compare our estimated luminosity function 
with data for satellit es in th e Loca l Group obtained from the 
sample compiled by Mateo (1998), supplemented with data 
from Irwin 2 . I Jnlike iMatecf l) 1998ft . we include the SMC, the 
LMC and M33 in our sample. Brighter than our estimated 
resolution limit of Mv < —12, we find that the luminosity 
function in our simulations is in excellent agreement with 
the Local Group data over a range of 7 magnitudes, down 
to Mv = —19. Of course, in the Local Group, the statistics 
at the bright end are rather poor: the last two data points 
(centred on —18 and —19.75) contain only one galaxy each 
(the LMC and M33, respectively). However, within the er- 
rors, the simulations are consistent with the data. 

We also compare o ur results in F i g.|2]wit h those in the 
semi-analytic model of iBenson et alJ i |2002rJL The simula- 
tions and the semi-analytic model are broadly in agreement 
over most of the luminosity range, from the resolution limit 
of the simulations to the brightest two bins. However, there 
is a significant difference at the bright end: while we find a 
satellite as bright as the LMC one third of the time, Ben- 
son et al. only find one such satellite 5% of the time. In 
the simulations, satellites as bright as M33 are produced in 
about 5% of systems, while the frequency in the Benson et 
al. model is less than 1%. The reasons for the disagreement 
between the semi-analytic model and the simulations (and 
Local Group data) are a reflection of the different treatment 
of various physical processes in the methods. However, the 
fact that our simulations match the observed satellite lumi- 
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Figure 2. The differential V-band luminosity function of satel- 
lites per central galaxy in our simulations. The blue error bars 
are lcr Poisson errors per magnitude bin. (For bins with more 
than one count the errors are simply vJV, but for bins with only 
one count, we take +viV as the upper errorbar and determine 
the lower errorbar by finding the mean of a Poisson distribu- 
tion whose integral from 1 to infinity corresponds to 16% of the 
distribution.) The Local Group satellite luminosity function, ob- 
tained from the compilation of Mateo (1998), supplemented with 
data from Irwin, is shown as red circles. The dashed line shows 
the median luminosity function of 70 realisations of Milky Way 
type halos calculated using a semi-analytic model by Benson et 
al (2002b) with the dotted lines indicating the 5-95 percentile 
width of their distribution. 



nosity function over a large range of magnitudes indicates 
that the relatively small observed number of satellites in the 
Local Group is not, in principle, difficult to explain within 
the CDM model. The much publicized "sate llite prob lem" 
exists, as other authors h ave remarked (e.g. iBullock et alJ 
120001 : [Benson et al.l2002hlL only when the astrophysical pro- 
cesses involved in the formation of visible galaxies are ne- 
glected. 
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Figure 3. Gas fraction in satellites as a function of V-band lumi- 
nosity. The satellites in the simulation are shown as squares and 
the brightest satellite in each system is surrounded by a red cir- 
cle. The Local Group satellites are shown as filled cyan circles and 
the LMC, SMC and M33 are indicated by crosses and labelled. 
Panel (a). Gas fraction as a function of the magnitude difference 
between the satellite and its central galaxy. For both simulations 
and data, only satellites with My < — 12 that lie within the virial 
radius of the central galaxy (the Milky way or M31 in the case 
of the real data) are plotted; in addition, for the real data a re- 
liable gas fraction measurement or upper limit is required. Panel 
(b). Gas fraction as a function of absolute V-band magnitude. 
Here, all Local Group members with reliable gas fraction mea- 
surements or upper limits are plotted. Local Group data are from 
Mateo (1998); SMC data are from Stanimorovi, Staveley-Smith 
& Jones (2003); LMC data are from Staveley-Smith et al. (2003); 
and M33 data are from McGaugh (2005). 



In order to investigate how realistic the satellites in our 
simulations are, particularly the bright ones, we compare, in 
Fig. |3] their gas fractions with those of real satellites. Panel 
(a) shows the gas fractions as a function of the magnitude 
difference between the satellite and the parent galaxy in the 
V-band. Only satellites with Mv < — 12 in the simulations, 
the Milky Way and M31 are shown and, for the latter, we 
require also a reliable measurement of, or upper limit to 
the gas fraction. The simulations show a trend of increasing 
gas fraction with decreasing luminosity which, as far as the 
scant data for the real satellites permit, seems consistent 
with the measurements. In particular, the brightest satel- 
lites in the simulations bracket the values measured for the 
LMC, the SMC, and M33 with 3 out of 9 simulated satel- 
lites having a larger gas fraction than the LMC and SMC. 
In panel (b) we extend the comparison by including not only 
satellites within the virial radius, but all the satellites in the 
Local Group that have reliable gas fraction measurements 
or upper limits. In this case, we plot the gas fraction against 
the V-band magnitude of the satellite. The trend of increas- 
ing gas fraction with decreasing luminosity is now clearer in 
the data and the locus they define agrees well with the lo- 
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Figure 4. Mean interior radial profiles for dark matter and satel- 
lites. The halo dark matter density profile, averaged over the 9 
systems with central galaxies, is shown by a red dashed line. The 
number density of satellites, also averaged over the 9 systems, is 
shown by a solid black line. The three blue lines show the num- 
ber density profiles of satellites for the three systems with 11 or 
more satellites (ghl, gh2 and gh3) and may be compared with the 
number density profiles for the 11 brightest satellites in the Milky 
Way (green squares and lines) and in M31 (magenta triangles and 
lines). In all these cases, the satellites were ordered according to 
distance from the centre and the density was evaluated at the mid 
point of the logarithm of the distance between adjacent satellites. 
The radial coordinate is in units of the virial radius, r v ; r , and all 
profiles have been normalised to the mean interior density at the 
virial radius. 

cus traced by the simulated satellite galaxies (although note 
that not all the points plotted are strictly speaking satel- 
lites, according to our definition, which requires satellites to 
lie within the virial radius.) We conclude that the satellites 
in the simulations not only have a luminosity function sim- 
ilar to that observed, but also have realistic gas fractions. 



4 THE SPATIAL DISTRIBUTION OF 
SATELLITE GALAXIES 

In Fig. 2] we show the mean interior radial mass density 
profile of the dark matter (red dashed line) and the mean 
interior radial number density profile of the satellites (thick 
black line) averaged over the 9 central galaxies in our simu- 
lations. Also plotted are the individual number density pro- 
files of the satellites in the three systems with 11 or more 
satellites (ghl, gh2 and gh3; blue lines), as well as the corre- 
sponding profiles for the 11 brightest satellites in the Milky 
Way (magenta line) and M31 (green line). All profiles have 
been normalized to their respective values at the virial ra- 
dius. The dark matter density pro file, shown as the dashed 
line, c losely follows the NFWform <|Navarro. Frenk fc White! 
Il996t iNavarro. Frenk fc White] Il997l) . As found in previ- 



Alignment of Satellite Galaxies 7 



ous N-body studies dStoehr et al l l2002t iGao et all 12004 
iLibeskind et al]l200fj) . the average satellite number density 
profile is flatter than the dark matter mass density profile. 
The three simulated systems with 11 or more satellites dis- 
play a range of profiles, but they all resemble the observed 
profiles of real satellites, and one of them (gh2), is remark- 
ably close to the Milky Way data over the entire radial range, 
0.1-lr vir . 

To characterize our simulations further, we investigate 
the flattening of the distributions of dark matter, satellites 
and stars in the central galaxy. We define the tensor of sec- 
ond moments, 

/< <■>■■:■ (i) 

where x^ is the j coordinate of the fiih particle in a ref- 
erence frame centred on the centre of mass of the central 
galaxy. To determine the flattening of the central galaxy, 
we consider only the innermost 98% of the stars and ex- 
clude the outermost 2% that make up the diffuse stellar 
halo. (Our results are insensitive to the exact fraction of ex- 
cluded stars, provided this is small.) To determine the prin- 
cipal axes of each distribution, we diagonalise the tensor 
Ijk- Its eigenvalues, a 2 ,b 2 ,c 2 , give the mean square devia- 
tion of the x, y and z coordinates along the principal axes. 
We define a > b > c as the major, intermediate and minor 
axes of the system, and use the notation {adm, 6dm, Cdm,}, 
{a sa t, b sat ; Csat }, {o g ai,6gal,£%al} to refer to the distributions 
of dark matter, satellites and stars in the central galaxy re- 
spectively. 

It is important to keep in mind that estimates of a, 
b and c using small numbers of objects are generally biased 
towards higher anisotropy. For example, if a system has only 
three members (e.g. satellite galaxies), our method would 
always return c = since three objects will always lie on a 
plane. Tests of the statistical robustness of our results are 
performed below. 

The Milky Way has 11 satellite galaxies within its virial 
radius with reliably measured positions and magnitudes. As 
shown in Table 1, of our sample of nine simulated central 
galaxies, three have 11 or more satellites and can be com- 
pared with the Milky Way. The axial ratios of each central 
galaxy in our simulations, its dark matter halo and, for sys- 
tems with 11 or more members within the virial radius, its 
satellite distribution are shown in Fig.|S] 

It is immediately apparent that, in the three sufficiently 
rich systems, the distribution of the 11 most massive satel- 
lites is considerably flatter than their respective dark matter 
halo, particularly in systems 1 and 3. Th i s is consistent with 
previous studies <|Libeskind et al.l 120051 iKang et alJ 120051 
IZentner et al.ll200ot) which also found large flattenings in 
simulated satellite systems, although in those studies the 
satellites were identified in dark matter simulations using 
a semi-analytic galaxy formation model rather than the full 
hydrodynamic calculations that we are analyzing here. Libe- 
skind et al ascribed this anisotropy to the preferential infall 
of substructures along the spines of filaments as they col- 
lapse to make a galaxy. 

In order to assess the statistical robustness of our re- 
sults, we performed three tests whose results are displayed 
in Table |1| For the first test, we constructed a spherically 
symmetric NFW halo with 10 6 particles. We then selected 
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Figure 5. Minor-to-major (c/a) versus intermediate-to-major 
(b/a) axial ratios for all the galactic systems in our simulations. 
Since a > b > c, the upper left triangular half of this plot cannot 
contain any points. Prolate objects lie along the diagonal while 
oblate objects lie along the right vertical axis. The numbers inside 
each symbol identify the galaxy. Axial ratios for stellar discs and 
for dark matter halos are shown for all nine systems by triangles 
and black squares respectively. For the three systems with 11 or 
more members, the axial ratios of the structure defined by the 11 
most massive satellites are shown by the circles. 



Table 2. The probability of randomly drawing values of the 
minor-to-major (c/a) and intermediate-to-major (b/a) axis ra- 
tio which are more extreme than those measured for the satellite 
systems in the simulations according to the three test described 
in the text. The three halos with 11 or more satellites within the 
virial radius are labelled ghl, gh2 and gh3 respectively. 





ghl 


gh2 


gh3 


Isotropic NFW sphere 


0.1% 


47.6% 


5.4% 


Squashed NFW profile 


0.4% 


51.5% 


6.4% 


Simulated dark halo 


0.6% 


66.2% 


10.5% 



11 particles at random from the halo 1000 times, and cal- 
culated the cumulative distributions of minor-to-major and 
intermediate-to-major axis ratios. The probability of draw- 
ing more extreme values that those measured for the satel- 
lite distributions in the simulations is given in the first row 
of Table 01 for each of the three galaxies with 11 or more 
satellites. For the second test, we construct flattened NFW 
halos by squashing a sphere according to the values of ad m , 
6dm, and Cd m for each of the three simulated halos. We then 
performed the same test as before; the results are displayed 
in the second row in Table. H Finally, we selected 11 dark 
matter particles directly from each of the three halos (again 
1000 times) and performed the same test with results given 
in the third row of the table. 
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Figure 6. The three simulated galaxies with 11 or more satellites. Clockwise, from the top-left, the plots show galaxies ghl, gh2 and 
gh3 respectively. The central galaxy, the satellite system and the dark matter halo are represented by the yellow, blue and red ellipsoids 
respectively. While each ellipsoid has been flattened and oriented according to the measured parameters, their relative sizes have been 
arbitrarily chosen for clarity. Note that the long, medium and short axes of the central galaxy are always aligned with the x, y, and z 



With the exception of galaxy gh2, whose satellite sys- 
tem is the least flattened of the three (see Fig[KJl, the num- 
bers in Table |1] rule out with high confidence the null- 
hypothesis that the satellite systems in the simulations 
were randomly drawn from either a spherical or a squashed 
NFW halo or from the dark matter distribution of the 



halo itself. This result agr ees well with those o btained by 
iLibeskind et alJ (l200rJ) and lZentner et alJ <2005^ . 

The flattening and relative orientations of the disc, halo 
and satellite system of each of our three well populated sim- 
ulations are shown schematically in cross section in Fig. |5] 
While the relative flattening and orientation of each compo- 
nent has been rendered according to the measured values, 
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for clarity, the relative sizes have been chosen arbitrarily. 
The central galaxy is represented by a yellow ellipsoid em- 
bedded within two larger ellipsoids, a blue one representing 
the satellite system and a red one representing the dark halo. 
The poles (i.e. the direction of the c axis) are visible as the 
point at which the lines of "longitude" converge. The various 
systems display a variety of orientations which we quantify 
below. 

The central galaxy orientations displayed in Fig. HJwere 
calculated from the innermost 98% of stars. To refine our 
estimate of orientation and for other purposes, we have per- 
formed a dynamical bulge-disc decomposition for each of 
the nine central gala xies in our s i mulat ions, applying the 
method propos e d bv lAbadi et all i200&t) and used also by 
lOkamoto et, al.l l|2005h . The net angular momentum, rela- 
tive to the centre of mass of the galaxy, of the material in 
the inner 10/i -1 kpc is used to define a "z" axis. The an- 
gular momentum, J z , about this axis is computed for each 
star particle and compared to the angular momentum of a 
circular orbit with the same energy, J C (E). All stars with 
Jn,/J C {E) > 0.75 are assigned to a disc component whose 
orientation is taken to be the direction of the nett angular 
momentum of the disc stars. (In some cases, this method 
may incorrectly assign some bulge stars to the disk, but 
this is not important for our analysis.) Generally, the direc- 
tion defined by the disc angular momentum coincides with 
that of the short axis of the overall stellar distribution, ex- 
cept in two case (galaxies gh8 and gh9) which are bulge- 
dominated galaxies with very small discs whose orientation 
is ill-defined. In what follows, we take the orientation of the 
galaxy to be that of the disc except for gh8 and gh9 for 
which we take the direction of their minor axis. 

The angles between the different components of each 
of the three galaxies with 11 or more satellites are plotted 
in Fig. |7| The top panel shows that in two of the galax- 
ies, the disc is inclined about ~ 45° relative to the minor 
axis of the dark halo. Surprisingly, in the third galaxy, the 
disc is orthogonal to the minor axis of the halo. The middle 
panel shows that in two cases (including the one with the 
orthogonal disc), the satellite systems are, within 20°, per- 
pendicular to the galactic disc. This is also surprising but 
it is exactly the alignment seen in the Milky Way galaxy, 
where the satellites lie approximately along a great circle on 
the sk y whose pole is in the galactic plane (e.g. lKroupa et all 
2005). In the third system (gh3), the satellite system is al- 
most perfectly aligned with the galactic disc. Finally, the 
bottom panel shows that in two of the three cases, the long 
axis of the satellite system is well aligned with the long axis 
of the dar k halo. This is consis tent with the conclusions 
reached bv lLibeskind et all i2005Tl from their (discless) dark 
matter simulations. The lack of alignment in system gh2 is 
probably due to the fact that this has the least aspherical 
dark halo whose long axis is therefore poorly determined. 



5 ALIGNMENT OF THE ANGULAR 
MOMENTA 

In this section, we extend our study of alignments to include 
the possible correlations between the angular momenta of 
the dark matter halo, J^, the satellite galaxy population, 
Jg^, and the stars of the central galaxy, J gal . Unless oth- 
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Figure 7. Alignments between the long (L) and short (S) axes 
of the different components of the three galaxies with 11 or more 
satellites. Panel (a): The cosine of the angles between the axes 
of the dark halo and of the galaxy. Panel (b): The cosine of the 
angles between the axes of the galaxy and of the satellite dis- 
tribution. Panel (c): The cosine of the angles between the axes 
of the dark halo and of the satellite distribution. The right-hand 
coordinates are labelled with the value of the angle. The most 
interesting relations are highlighted with stars. 



erwise noted, we consider a number-weighted mean angular 
momentum, J = (rxv). For the dark matter halo and the 
central galaxy whose particles all have roughly the same 
mass, this "specific" angular momentum is approximately 
proportional to the standard, mass-weighted angular mo- 
mentum ((rxp)). For a satellite system, however, this def- 
inition gives equal weight to each satellite, preventing the 
statistics from becoming dominated by one or two very large 
satellites. 

The angle between j dm and J gal evaluated for all the 
material within the virial radius, is plotted in Fig. |H| Since 
the number of satellites is immaterial in this case, we expand 
our sample to include the six other galactic systems (identi- 
fied according to the criteria outlined in Sec. 12.21 that were 
excluded from the satellite analysis in the preceeding sec- 
tion. The angular momenta of the galaxy and the dark mat- 
ter within the virial radius are aligned to within ~ 30°, i.e. 
the galaxy spins in essentially the same direction as the dark 
matter halo. The single exception is system gh9 in which the 
galaxy and the dark matter are counter-rotating. A KS test 
shows that the probability of obtaining the distribution of 
cos(9 shown in the figure from nine objects drawn at random 
from a larger sample of randomly oriented galaxies is only 
~ 0.016%. We conclude that there is a significant alignment 
between the angular momenta of the galaxy and that of the 
dark halo it inhabits. This is consistent with the acquisition 
of angular momentum prior to the collapse of the system, 
when dark matter and gas we re well mixed , as expected in 
the tidal torque theory (e.g. IWhitd Fl984h and as usually 
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Figure 8. Alignment of the angular momentum vectors of the 
galaxy (J ga j) and the dark halo (J<j m ) for material within the 
virial radius for our nine simulated galaxy systems. The right- 
hand axis gives the value of the angle. Each galactic system has 
been assigned a colour that will be retained throughout the re- 
mainder of this paper. 



assumed in sem i-analytic models of galaxy formation (e.g. 
ICole et aljEiooT) . 

The radial dependence of the galaxy-halo spin align- 
ment is shown in Fig. |5] where we plot the cosine of the 
angle between J gal and J dm (< r) for our nine systems. In 7 
out of 9 cases, the alignment persists out to the virial radius 
of the halo, never straying much beyond ~ 30°. The two 
exceptions are systems gh9 and gh8. In the first of these, a 
tight alignment persists for r < 0.7r v i r , but it rapidly disap- 
pears beyond this radius. A visual inspection of this system 
shows that the reason for the rapid change in spin direction 
is simply the presence of a large counter-rotating fragment 
of dark matter in the outer parts of the halo which has re- 
cently been accreted, causing the halo angular momentum 
to flip in the outer parts. The other anomalous system owes 
its strange behaviour to the fact that the galaxy is almost 
a pure spheroid with very little angular momentum and for 
which the direction of J gal is poorly defined. 

Previous studies (e.g. IWarren et alj Il992l 
iBailin fc Steinmetzl I2005T) have found that the halo 
angular momentum, J^, tends to be alig ned with the 
halo minor axis, c dm . The simulations of IWarren et alJ 
(1992) clearly show such an alignment in the inner parts 
of the dark halo. Bailin & Steinmetz (2005; see their 
Fig. 16) demonstrate that the alignment becomes weaker 
with radius but is still significant (~ 30°) at the virial 
radius. By contrast, IBailin et alJ i2005l) report an absence 
of alignment between c dm and J gal at the virial radius 
in seven hydrodynamic simulations. We now consider the 
angle subtended by J dm and c^ in our own simulations. 

In Fig. llOl we plot the angle between J dm and c dm for all 
the material within the virial radius and, in Fig. 1111 the run 



Figure 9. The angle between the angular momentum of the disc 
GLjisc) an< ^ ^he dark matter halo (J dm (< r)) as a function of 
radius. Each colour represents a single system according to the 
same code used in FiglSl 
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Figure 10. The angle between the angular momentum and the 
short axis of the dark matter halo (J dm . c dm , pentagons) and 
of the galaxies (J ga j • £ ga ii squares) for the nine galactic systems. 
The angular momentum vector in most of the galaxies points in 
the direction of the shortest galactic axis. For the halos, however, 
no such correlation exists. 




Figure 11. The angle between iLj m and Cd m for material within 
a given radius. Each of the 9 galactic halos is plotted using a 
different colour, according to the same code used in Fig. 8. 

of this angle with radius. We also plot the KS probability 
that the distribution of angles in the nine systems at each 
radius is consistent with a uniform distributio n. Our results 
are qu alitat ively in good agreem e nt wit h both lWarren et ail 
ljl992h and iBailin fc Steinmet,3 l|2005h . The radial depen- 
dence of J dm . c dm , displayed in Fig. ^3 shows that, apart 
from one anomalous system, gh9, there is a good alignment 
between these two vectors in the inner parts of the halo, 
within a few tenths of r v i r . In the outer parts the correlation 
becomes much weaker, which is perhaps not surprising for 
systems whose shape is supported by an anisot ropic veloc- 
ity d i spersion tensor rath er than by rotation jFren k et alJ 

[ 19851: IWarren et alJll992h While IWarren et alJ <1992h and 
Bajhr^^BtelmneM 12005T) find a weak correlation at the 
virial radius our small sample of 9 halos is consistent with 
no correlation. Thus, the KS probability that the cosine of 
the angle for these 9 halos is consistent with a uniform dis- 
tribution jumps from ~ 15% at > 0.7r v i r to ~ 80% at r v i r . 

Turning to the galaxies, Fig. 1101 shows that in all but 
two of them, the angular momentum vector points in the di- 
rection of the shortest galactic axis, as expected for systems 
flattened by rotation. Of the two discrepant galaxies, gh8 is 
a nearly perfect spheroid with very little rotation while gh9 
is a rapidly rotating spheroid, whose angular momentum is 
dominated by the bulge's bulk rotation. 

The orientation of the galaxy and its host halo is dis- 
played in Fig. ^] where we plot the angle between c dm and 
j gal as a function of radius. Since J gal is always parallel to 
fidisc, this is equivalent to plotting the angle between c dm 
and Cdisc- Apart, again, from the anomalous gh9, the angu- 
lar momentum vector of the galaxy points along the short 
axis of the halo in the inner parts of the system. This align- 
ment begins to weaken beyond ~ 0.5r v i r , but even at r v i r 
the KS probability is consistent with a uniform distribution 
only at the ~ 15% level. Our results agree well with those 
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Figure 12. The angle between the angular momentum of the disc 
(Jgal) and the short axis of the dark matter halo (c<i m ) for the 
9 galactic system as a function of radius. Each of the 9 galactic 
halos is plotted using a different colour, according to the same 
code used in Fig. 8. 

of IBailin et alJ i2005l) in the inner parts of the halo, but 
they are marginally inconsistent at the virial radius where 
their small sample of objects has a distribution of cosines 
consistent with uniform. 

We now turn to the alignment of the angular momen- 
tum of the satellite system with those of the galaxy and the 
dark matter halo. In this case, we can use only the three 
systems with a large satellite population. The results are 
plotted Fig. I13f a) as filled and empty squares. For all three 
systems, there is some alignment (within ~ 40°) between 
the spin vectors of the satellite systems and the dark halo. 
However, in only one case (gh3) is the spin of the satellites 
aligned with that of the galaxy. 

Finally, we ask whether the net angular momentum of 
the satellite system is related to the shape of the system, 
i.e. whether the satellites tend to orbit in the approximate 
plane that they define. If this were the case, we would ex- 
pect j aat to point along the short axis of their distribution, 
£saf I n f ac t, this is not what we find. J sat tends to point in 
a direction which, while often approximately at right angles 
to a sat , the long axis of the satellite system (Fig. I13f b)). 
it is nevertheless not along the short axis of the system, 
c sat . This arrangement is illustrated schematically in Fig. 1141 
where the blue ellipsoid demarks the structure defined by the 
anisotropic satellite distribution, while the spin axis, , is 
shown as a red line and its three projections along the three 
principal axes of the satellite structure are shown as black 
lines. From Fig. Il^f b) we find that in all three cases the 
angular momentum of the satellite system subtends an an- 
gle greater than ~ 60° with the long axis of the satellite 
distribution, indicating that the satellites are not orbiting 
in the apparent plane defined by their spatial distribution. 
However, there is little consistency amongst the 3 galactic 
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Figure 14. Schematic diagram showing the satellite system angular momentum vector, J~ at , in three dimension (red line) calculated as 
(rxv). The blue ellipsoids illustrate the shape of the satellite distribution. The black lines are the projections of the vector J gat along 
the three principal axes of the satellite distribution. Clockwise, from the left, we show galactic systems ghl, gh2 and gh3 
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Figure 13. Top panel (a): the cosine of the angle between the 
satellite angular momentum vector, J sat , and J dm (filled squares) 
and between and J gal (empty squares). Bottom panel (b): 
the cosine of the angle between J sat and the long axes of its 
spatial distribution, |L at . The filled squares correspond to the 
case when the angular momentum is calculated assuming each 
satellite has the same mass (i.e J = (r X v)), while the empty 
squares correspond to the case when the angular momentum is 
calculated in the standard way (i.e J = (r X p)). 



systems which exhibit a variety of orientations and align- 
ments. 



6 CONCLUSION AND DISCUSSION 

We have analyzed two N-body/SPH simulations of galaxy 
formation in a ACDM universe, one of a single bright galaxy 
and the other of a small region. Our simulations include 
the main physical processes thought to be important in 
galaxy formation: metal-dependent gas cooling, heating by 
photoionizing the primordial hydrogen gas, star formation, 
metal production and feedback due to supernovae energy 
injection. In total, we obtained a sample of 9 bright galax- 
ies comparable in luminosity to the Milky Way, containing a 
total of 78 satellites. Three of the central galaxies had a pop- 
ulation of at least 11 satellites, the number of well-studied 
satellites in the Milky Way. We have used these samples 
to investigate several interrelated properties of the satellite 
population. 

We first investigated the luminosity function of satel- 
lites, a property often r egarded as a challe nge to the cold 
dark matter model (e.g. iMoore et al"Hl9 99). We find, how- 
ever, that our simulations produce an excellent match to 
the observed luminosity function of satellites in the Local 
Group, at least to the resolution limit of the calculation 
which corresponds to My- ~ —12, close to the luminosity 
of the faintest observed satellites. The gas fractions in these 
satellites agree well with observations, suggesting that their 
formation paths may be similar to those of real satellites. 
Many of the small subhalos resolved in the si mulations fail 
to mak e a substantially bright satellite. As iBenson et alJ 
(2002b) have argued, star formation is inefficient in small 
subhalos due to the combined effects of reionization and su- 
pernovae feedback which limit the supply of cool gas to the 
subhalo. 

The match of our simulations to the satellite data over 
most of the observed range suggests that the relative paucity 
of satellites in the Local Group (the "satellite problem") 
does not reflect the much larger abundance of subhalos, but 
rather feedback e ffects that limit the gro wth o f small g alax- 
ies. F urthermore, IBenson et al] <l2002bl) and IStoehr et al] 
( 2002) have argued that the internal dynamics of CDM satel- 
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lites are also consistent with the data. Overa ll, our results 
agree w ell with the semi-analytic model of iBenson et alJ 
(2002b). In the semi-analytic model, however, bright satel- 
lites, of luminosity similar to the LMC or M33, were rare. 
This is not the case in our simulations which match the 
bright end of the satellite luminosity function well. The dif- 
ferences between the two treatments are no doubt due to the 
different ways in which the physics of galaxy formation are 
modelled, including differences in the treatment of feedback 
and star formation and in the satellite merger rates. 

We next considered the spatial distribution of the satel- 
lites. In the outer half of the system, the number density pro- 
file of the satellites tracks the dark matter well but in the 
inner parts it falls below the dark matter profile. Remark- 
ably, the number density profile of the 11 brightest satellites 
in one of our simulations is almost indistinguishable from the 
profile for the 11 brightest satellites of the Milky Way. The 
profiles of the other 2 well populated systems are broadly 
similar, as is the profile fo r the satellites of M31. As in pre- 
vious N-body studies ( e.g. iGao et al"ll2004l : iLibeskind et alJ 
l200fJ; ISh aw et af]l2006T) . we find no resolved satellites in the 
inner ~ 10% of the halo. 

We can study the distribution of individual satellite 
systems, particularly their alignments with the halo and 
disc, for the three systems that formed more than 11 satel- 
lites. In agreement with previous purely N-body studies 
iLibeskind et all2005llZentner et a.lj|2005ll . we find that the 
satellites tend to be distributed in a highly flattened config- 
uration whose major axis is alig ned with the major ax is of 
the (generally triaxial) dark halo. ILibeskind et al.l J2005i) ar- 
gued that this arrangement reflects the preferential infall of 
satellites along the spine of the filaments of the cosmic web. 
Our gasdynamic simulations allow us to go further than pre- 
vious N-body work and investigate the alignment of satellite 
systems with galactic discs and the alignment of the latter 
with the dark matter halo. We find that in 2 out or 3 sys- 
tems, the satellite system is nearly perpendicular (to within 
20°) to the plane of the galactic disc. This surprising config- 
uration is exactly what is seen in the Milky Way. The third 
satellite system ended up well aligned with the disc. 

To investigate these alignments further, we calculated 
the relationship between the disc and the halo, this time 
using the 9 bright galaxies in our simulations. Previous hy- 
drodynamic simulations have found a good alignment be- 
tween the disc and the principal plane of the halo (i.e. 
a good correlation between the directions of the disc axis 
and the hal o minor axis), but only in the inner ~ .2rvi r 
of the halo (iKazantzidis et aljl20fj3 : iBailin et al.ll2005l) . Be- 
yond this, these studies find little or no correlation. We also 
find a good alignment in our simulations between the disc 
and the halo at small radii but, unlike in the previous stud- 
ies, the correlation persists, albeit much weakened, out to 
the virial radius. This sort of alignment is perhaps not un- 
expected in the simplest interpretation of the tidal torque 
theory since both the dark halo and the baryon component 
experience similar tidal torques. The re latively small dif- 
ferenc es be tween our results a nd those of IKazantzidis et alJ 
( 2004) and lBailin et alJ l|2005fi are most likely the result of 
different treatments of star formation and feedback. For ex- 
ample, the strong feedback in our simulations at early times 
which leads later to a prolonged period of gas accretion and 
the formation of a large disc in a relatively quiet halo favours 



the persistence of a relationship between the properties of 
the halo and the disc. 

Of the two cases in which the satellite systems are nearly 
perpendicular to the disc, in the one which exhibits the 
strongest alignment (ghl) the disc is also nearly perpen- 
dicular to the minor axis of the halo. In the other case, 
the alignment between the disc and the halo is weaker. In 
the third system (gh3), in which the satellite system lies in 
the plane of the disc, the disc is only roughly aligned with 
the halo. Although our sample is small, it suggests that it 
should not be surprising that similar kinds of alignments are 
found in observational studies of satellites in the SPSS and 
2dFG RS (lBrainerj200llYang et ai1l2005l ISaTes fc Lambad 
l2004h . 

Finally, we investigated the connection between the an- 
gular momenta of the disc and the dark halo. For most sys- 
tems, the two vectors are very well aligned as a function of 
radius, out to the edge of the system. The halo spin tends 
to be point along the short axis of the halo. However, since 
the shape of the halo varies with redshift, the spin axis of 
the disc or of the halo can often loose its alignment with 
the short halo axis. Thus, care should be taken when in- 
terpreting alignment statistics, either in simulations or in 
observational data. 

In conclusion, our simulations have revealed a number 
of interesting connections between the properties of central 
galaxies, their satellite systems and their dark matter halos. 
They also point to a good overall match between predic- 
tions of the ACDM model and the observations of galaxy 
satellite systems that we have considered here such as the 
luminosity function and the peculiar spatial arrangement of 
the satellites of the Milky Way. However, due to their high 
computational cost, our simulations are still too small to 
provide good statistics. Larger simulations of this kind are 
required to test the validity of the trends that we have found. 
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